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Abstract 

We review the extension of the factorization formahsm of perturbative QCD to soft 
rescattering associated with hard processes in nuclei. 

1 Introduction 

This paper reviews a perturbative QCD (pQCD) treatment of the hard scattering of hadrons 
and leptons in nuclei, based on factorization, It describes in part work in collaboration with 
Ma Luo []T], 1^, ^ and Xiaofeng Guo [|, §,0. At the outset, it may be useful to clarify 
the relation of this work to the works of Baier et al. (BDMPS) [0 and Zakharov P). For 
more recent progress in relating these two approaches, see §]. We have tried to illustrate 
this relation schematically in Fig. |l|. The BDMPS analysis begins (Fig. |l]a) with the classic 
treatment of radiation induced when a charged particle passes through a large target, due 
originally to Landau, Pomeranchuk and Migdal (LPM). This analysis does not require the 
presence of a hard scattering, but describes the coherent results of many soft scatterings. 
Its primary subject has traditionally been induced energy loss. Our analysis (GLQS) begins 
with the perturbative QCD treatment of hard-scattering in a relatively small target (Fig. 0b), 
in which the primary subject of interest is momentum transfer. A complete analysis (Fig. |I|c) 
of hard scattering in a large target involves both energy loss and the transverse momenta due 
to initial- and final-state soft scatterings. Our work is a step in this direction, attempting 
to stay as close as possible to the pQCD formalism, in which we may readily quantify 
corrections. To be specific, we consider only a single soft initial- or final-state interaction in 
addition to the hard scattering. Our central observation is that for suitably-defined jet and 
related inclusive cross sections this is the first order in an expansion in the quantity 

where A represents a nonperturbative scale, which we shall identify with a higher-twist parton 
distribution below. That additional scatterings are suppressed by factors of is perhaps 
surprising. Let us review why this is the case, at least for certain cross sections. 
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Figure 1: Alternate approaches to hard scattering in nuclei, (a) Landau-Pomeranchuk- 
Migdal analysis treats energy loss due to many soft scatterings, (b) Perturbative QCD 
analysis treats momentum transfer due to hard scattering, (c) For scattering in nuclei, both 
must be combined. 

The basic analysis of hard-scattering in nuclear matter (cold or hot) |T0[ is quite simple. 
To be specific, consider the scattering of a quark, as shown in Fig. ^. A hard-scattering with 
momentum transfer Q can resolve states whose lifetimes are as short as for instance 
quarks off-shell by order Q, but still less that Q. The off-shellness of the scattered quark 
increases with the momentum transfer simply because the number of available states increases 
with increasing momentum. Similarly, the scattered quark, of momentum p' is typically off- 
shell by order nij < Q. We may think of mj as the invariant mass of the jet into which 
quark fragments. If we are to recognize the jet, we must have mj -C Ej = p'^, with Ej 
being energy of the jet. On the other hand, the counting of available states ensures that 
mj > Aqcd- 

Now the scattered quark has a lifetime in its own rest frame At^^ ^ ~ ^ with mj -C 
Ej. In the target rest frame, however, this becomes, for large enough Ej/mj, At^^^^^'^^^ ~ 
m7 (mj) ^ "^here Ra is the (fixed) target size. Thus, at high enough energy the lifetime 
of the scattered quark will exceed the target size, even though the quark itself is far off the 
mass shell, typically by a scale that grows with the momentum transfer Q. 

Further couplings of the off-shell quark are suppressed, first of all by the strong coupling 
evaluated at scale mj, and, more importantly, by an overall factor of ~ l/Q^, since 

the effective size of the scattered quark decreases with momentum transfer in this manner. 

In summary, for inclusive processes such as jet production, high-Q implies that process- 
dependent multiple scattering is power-suppressed compared to single scattering. The most 
important point here is that the scattered particle remains off-shell for its entire transit 
of the target. Thus, its interactions with the target may be treated by the formalism of 
perturbative QCD, which, however, must be extended to include corrections that decrease 
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Figure 2: Sketch for the scattering of 
a quark of momentum xp in a large 
nucleus. 



Figure 3: Classification of parton multiple scatter- 
ing in nuclear medium: (a) interactions internal to 
the nucleus, (b) initial-state interactions, and (c) 
final-state interactions. 



with extra powers of momentum transfer. Up to the first such "higher-twist" contribution, 
a general cross section has the representation [|ll] 

a(g) = if°®/2®/2+ (^^)i/^®/2®/4 + 0(^^^ , (2) 

where ® represents convolutions in fractional momenta carried by partons, and represents 
a parton distribution of twist n. Target-size dependence due to multiple scattering can only 
appear in the second term in this expansion. 



2 Parton-Nucleus Scattering in Perturbative QCD 

To distinguish parton-nucleus multiple scattering from partonic dynamics internal to the 
nucleus, we classify the multiple scattering internal to a nuclear target in the following 
three categories: (a) initial-state interactions internal to the nucleus, (b) initial-state parton- 
nucleus interactions (ISl), and (c) final-state parton-nucleus interactions (FSI), as shown in 
Fig. |[ To a certain degree, this classification is ambiguous, but it can be made well-defined 
if we are careful. 

Initial-state interactions internal to the nucleus change the twist-2 parton distributions 
of the nucleus, as shown in Fig. |^a. As a result, the effective parton distributions of a large 
nucleus are different from a simple sum of individual nucleon's parton distributions. This 
is known as the "EMC" effect for the region where the parton's momentum fraction x is 
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Figure 4: Scattering amplitude between 
two incoming hadrons with one hard col- 
hsion. 



Figure 5: Perturbative QCD factorization 
at leading twist. 



not too small. Since only a single parton from the nucleus participates the hard collision to 
leading power, the effect of the initial-state interactions internal to the nucleus is (almost by 
definition) leading twist. The A-dependence of the "EMC" effect provides a relatively small 
nuclear size dependence to the first term in the Eq. @. 

On the other hand, the initial-state and final-state parton-nucleus interactions, as shown 
in Fig. and |^c, involve at least two physical partons from the nucleus at the hard collisions. 
Thus the "Cronin effect", ^"-dependence with a > 1, due to multiple scattering, is higher- 
twist for inclusive distributions. 



2.1 Factorization at Leading Powers 

Let us review some of the details of a factorized cross section like the one in Eq. (^. 

The leading power contributions to a general cross section involve only one hard collision 
between two partons from two incoming hadrons (or nuclei), as shown in Fig. ^ An energetic 
third parton is produced in the collision, which fragments into either a jet or a hadron of 
momentum p'. After squaring the scattering amplitude in Fig. |^, and summing over all 
final-states, the cross section can be factorized into a form like the first term in Eq. (H), 




{xipi,x'p2,p'/z) , (3) 



where J2ijk ^^ns over all parton flavors and all scale dependence is implicit. The D^/k are 
fragmentation functions for a parton of type k to produce a hadron h. For jet production, the 
fragmentation from a parton to a jet, suitably defined, is calculable in perturbation theory, 
and may be absorbed into the "hard scattering function" a. Then, the factorized single 
scattering formula in Eq. (|) is reduced to the first term in Eq. (Q). The fa/p are twist-2 
distributions of parton type a in hadron p. They have the interpretation of expectation 
values in the hadronic state of products of fields on the light cone, for instance, for a quark 
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distribution 

Mx,Q) = J ^e-P"^-(p|g(0)^g(y-)b) , (4) 

where for simplicity we choose the A'^ = gauge, assuming p is in the plus direction. Eq. (j^) 
is illustrated by Fig. ^ As shown, the convolution in Eq. (j^) is in terms of the momentum 
fractions Xi and x' carried by partons i and j, from hadrons pi and p2, respectively, into the 
hard scattering. 

The factorized formula in Eq. (|^) illustrates the general leading power collinear factoriza- 
tion theorem |]T^. It consistently separates perturbatively calculable short-distance physics 
into a, and isolates long-distance effects into universal nonperturbative matrix elements (or 
distributions), such as fa/p or Dh/k, associated with each observed hadron. Quantum inter- 
ference between long- and short-distance physics is power-suppressed, by the large energy 
exchange of the collisions. Predictions of pQCD follow when processes with different hard 
scatterings but the same nonperturbative matrix elements are compared. 

In the case of collisions on a nuclear target, the factorized single scattering formula 
remains valid, except that the twist-2 parton distribution fa/p is defined on a nuclear state, 
instead of a free nucleon state. For example, for a nucleus of momentum Pa = Ap^, the 
effective quark distribution is defined in the same way as in Eq. (^), with \p) a nuclear state, 
\Pa)- Such an effective nuclear parton distribution includes the "EMC" effect, and is still a 
twist-2 distribution function by the definition of its operator. 

Power-suppressed corrections to Eq. involve ratios of the nonperturbative momentum 
scales in the hadron, A ~ Aqcd ~ 1/fm, over the energy exchange of hard collisions, Q, as 
(A^/Q^)". These corrections can come from several different sources, including the effect of 
partons' non-coUinear momentum components, and the effects of interactions involving more 
than one partons from each hadron. In the case of nuclear collisions, such power-suppressed 
collisions can be enhanced by the nuclear size, as in Eq. (p. 



2.2 Factorization at Nonleading Powers 

Much of the predictive content of pQCD is contained in factorization theorems like Eq. (|^). In 
order to consistently treat the power suppressed multiple scattering, we need a corresponding 
factorization theorem for higher-twist (i.e., power suppressed) contributions to hadronic hard 
scattering. 

Fig. ^ is a picture for a power suppressed contribution to hard scattering. In this case, 
two partons i and i' with momenta Xipi and X2P1 from the target (the "nucleus" ) collide with 
a single parton j of momentum x'p2 (from the "projectile"). After squaring the scattering 
amplitude in Fig. |, and summing over all final states, the power suppressed contribution to 
the cross section can be factorized into the form [|TT|] 

E J^D,/,iz) (5) 

which can be illustrated by the sketch in Fig. |^. The expectation value T corresponding to 
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Figure 6: Scattering amplitude between 
two incoming hadrons with two hard coUi- 
sion. 



Figure 7: Perturbative QCD factorization 
at the next-to-leading power. 



this multiparton contribution from the target is typically of the form 



T^ii')/p{xi,X2,X3,Q) 



(7) 

where Bi is the field corresponding to a parton of type i = q,q,G. In Eq. (^), the hard 
part depends on the identities and momentum fractions of the incoming partons, 

but is otherwise independent of the structure - in particular the size - of the target (and 
projectile). To find A-enhancement due to multiple scattering, we must look elsewhere. 

Before identifying the source of the A-enhancement, we briefly explain why the factor- 
ization formula in Eq. can be valid. Although the formal arguments for the validity 
are well-documented [[Tl|, a heuristic understanding of factorization can be useful. Such a 
understanding may be found in the Lorentz transformation properties of gauge fields |]T3|. 

In hadron-hadron collisions, the factorization could be broken if interactions of long- 
range fields, labeled by "5"' in Fig. H, between the two incoming hadrons are important. The 
interactions of the long-range fields could alter the hadronic states of the incoming hadrons, 
and subsequently, change the parton matrix elements (or distributions) before the hard 
collisions take place. Without universality of the parton matrix elements (or distributions) 
in Eqs. (|) and (H), the factorized formulas would lose predictive power, and we would say 
that factorization is broken. Formal proofs of the factorization theorem must show that all 
such long-range soft interactions are either power suppressed, or can be effectively removed 
due to unitarity, for inclusive observables [|r^. 

Let us now review a heuristic argument that soft interactions between the two incoming 
hadrons are kinematically suppressed, due to the Lorentz transformation properties of gauge 
fields [0. Consider a collision between hadron A and hadron B, as shown in Fig. |^. Hadron 
B is moving very fast along the direction of 0:3 with /3 ~ 1 , while hadron A is at rest in the 
x' frame. Let us define A = f3ct' — x'^, where is the third component in the x' frame. 
The long-range fields generated by hadron B, at rest in the x frame, behave very differently 
in the x and x' frames, and different types of fields have different properties under Lorentz 
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Figure 8: Interactions of long- 
range fields that might break the 
factorization. 



Figure 9: Sketch for soft interactions of long-range 
fields between two colliding hadrons A and B. 



transformations. For example, the three types of fields: scalar field, longitudinal component 
of a gauge field, and its E3 field strength, have the following behaviors; with 7 = (1 — /?^)~^/^, 
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(9) 



(Xf. + 72A2)3/2 



=^ — "strongly contracted!" (10) 

T 

Although the magnitude of the longitudinal component A'^{x') of the gauge field is not 
suppressed under the Lorentz transformation, as shown in Eq. (^, and its interactions can 
be very strong, a short calculation shows that as /? — *• 1, it becomes gauge equivalent to a 
vanishing gauge field. On the other hand, the gluon field strength is suppressed under the 
Lorentz transformation even more strongly than the scalar field, Eq. (^. In terms of energy 
scales, the 1/7^ in Eq. (|T0|) translates into a suppression factor of l/Q"^, which suggests 
that the factorization should be valid at the order of 1/Q^, and might fail at 0(1/ Q"^) 



11, 13, 14, 15 



Showing the factorization at the next-to-leading power is a beginning toward a unified 
discussion of 0{1/Q^) effects in a wide class of processes. A systematic treatment of double 
scattering in a nuclear medium is an immediate application of the generalized factorization 
theorem. 
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Figure 10: Poles that give rise to an A-enhanced cross section in deeply inelastic scattering. 
2.3 ^4- Enhancement from Matrix Elements 

As we pointed out in last subsection, the partonic hard part in the factorized formula in 
Eq. (|^) is independent of the structure - in particular the size - of the target. Therefore, 
we need to find the A-enhancement due to multiple scattering from the matrix elements, if 
there is to be any size enhancement. 

For definiteness, we consider photoproduction or deeply inelastic scattering on a nucleus 
fl], ^. In this case, the additional soft scattering is always a final-state interaction. The 
structure of the target is manifest only in the matrix element T in Eq. (^. Each pair of fields 
in the matrix element Eq. (0) represents a parton that participates in the hard scattering. 
The y~ integrals parameterize the distance between the positions of these particles along the 
path of the outgoing scattered quark. In Eq. (^, integrals over the distances y~ generally 
cannot grow with the size of the target because of oscillations from the exponential factors 



e 



Since the kinematics of a single-scale hard collision is only sensitive to the total momen- 
tum from the target, two of the three momentum fractions: Xi,X2, and X3 cannot be fixed 
by the hard collisions. If the integration of the momentum fractions is dominated by the 
region where X2 ~ and X3 ~ 0, the corresponding y~ integration in Eq. (0), 



j dy^ e 



oc size of the target. 



is proportional to the size of the target or the A^^^. 

The question now is if it is possible that the partonic parts, (y^^l)j^j in Eq. (||) can have 
contributions that are dominated by regions where one or two of the three momentum 
fractions vanish. The answer is yes. Some of the ^"[jy)^^ in the Xi integrals have two poles, 
labeled by the crosses "x" in Fig. ID. The symbols "x" in Fig. ^ represent the potential 
poles when the corresponding propagators become on-shell. In the example of deeply inelastic 
scattering, the poles are always associated with the scattered particle |^, while in other 
processes, the poles can also be associated with incoming particles 0, |[. It is important 
to emphasize that using a pole in the complex Xi (longitudinal momentum) space to do 
the integral does not correspond to assuming on-shell propagation for the scattered quark. 
Indeed, the Xi integrals are not pinched between coalescing singularities at such points, and 
the same results could be derived by performing the Xj integrals without ever going through 
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Figure 11: Sketch of the twist-4 quark- 
gluon correlation function that gives rise 
to an A^^^ type enhancement. 



Figure 12: Sketch of the twist-4 gluon- 
gluon correlation function that gives rise 
to an A^^^ type enhancement. 



the Xi = points. It is also worth noting that this is a feature unique to higher-twist matrix 
elements, for which the Xi are not restricted to be positive definite. Physically, this is possible 
because, unlike leading twist matrix elements, they do not generally have the interpretations 
of probabilities. 

The result of this reasoning is that two of the three momentum fraction integrations: 
dxi dx2 dx3 in Eq. (|^) are fixed by the two poles, and the convolution over dxi dx2 dx^ in 
Eq. (Bl) is simplified to an integration over only one momentum fraction. 



dxi dx2 dx3 T(^ii')/pi (xi, X2, x^) a^f-,)^j{xipi, x'p2,p'/ z) 
^ j dxTq{x,A) a^^^^-{xpi,x'p2,p' /z), 



(11) 



where the partonic part a(^) is finite and perturbative with the superscript [D) indicating 
the contribution from double scattering. The above matrix element, Tq{x,A), as illustrated 



in Fig. [TT|, has the form 

T,{x,A) 



27r J 1-n 

X ^ kVA\mi^F^^{y^)F^Ay2)(l{Vi)\VA) 



(12) 



where \pa} is the relevant nuclear state. The variable x here is the fractional momentum 
associated with the hard parton from the target that initiates the process. The soft scattering 
contributes a negligible longitudinal fractional momentum. Details of the reasoning and 
calculation for deeply inelastic scattering are given in Ref. [Q. 

Similar to the quark-gluon correlation function Tg(x, A) in Eq. (|T2|), another important 
twist-4 parton correlation function that gives rise to an A}!"^ enhancement is the gluon-gluon 
correlation function [lltl, 



T,{x,A) 



27r 



dy2 di/3 
27r 



X 



xp^ 



{pA\F^nO)F-^{y,)F\{y,)F\iy^)\p^) 



(13) 



which is illustrated in Fig. 112 
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In this form of the twist-4 parton-parton correlation functions, two integrals over the 
y~ and can grow with the nuclear radius as fast as A^/'^. However, if we require local 
color confinement, the difference between the light-cone coordinates of the two field strengths 
should be limited to the nucleon size. Therefore, only one of the two y~ integrals can be 
extended to the size of nuclear target. The twist-4 parton-parton correlation functions are 
then proportional to the size of the target, that is, enhanced by A^^^. 



3 Applications 

In Refs. |]1| and , we have applied the formalism sketched above to single-particle inclusive 
and single-jet production for deeply inelastic scattering and photoproduction. These cases 
involve final-state interactions only. In each case, the leading l/Q"^ correction is proportional 
to the matrix elements in Eqs. ([T2|) and (pIS]). Of course, the value of the correction cannot 
be estimated without an idea of the magnitudes of the T's. Since these magnitudes are 
nonperturbative they must be taken from experiment. At the same time, we expect the 
x-dependence of the probability to detect the hard parton to be essentially unaffected by 
the presence or absence of an additional soft scattering. Thus, we choose the ansatz |l[ 

T,(x,A)=AMi/3/,/p^(a;,A) (14) 

for i = q,g in terms of the corresponding twist-two effective nuclear parton distribution fi/^, 
with A a constant with dimensions of mass (see Eq. (|l])). This assumption facilitates the 
comparison to data. 



3.1 Momentum Imbalance of Di-Jets in Photoproduction 

A quantity that is sensitive to final-state rescattering in a particularly direct way is the 
momentum imbalance of di-jets in photoproduction in nuclei. The A^/^ dependence of this 
quantity measured by the Fermilab E683 collaboration is a clear signal of the presence of 
double scattering |]16| . 

The di-jet momentum imbalance kT^p measured at Fermilab is defined as [0 



= ^ (pti + PT2) sin A0 , (15) 
where and are the transverse momenta of the di-jets and the angle between the di-jets 



A0 is defined in Fig. |I3|. This momentum imbalance is the projection of the momentum 
imbalance out of the scattering plane defined by the beam axes and either of the two observed 
jets. The data from E683 shows that the averaged momentum imbalance (k^^) can be 
parameterized as (1.44 -|- 0.174^4^/^) GeV^. The term proportional to A^^^ should be a 
consequence of the multiple scattering of these two jets when they pass through the nuclear 
medium 0. 

In Ref. 0, we calculated the di-jet momentum imbalance in pQCD by applying the 
twist-4 factorization formalism to double scattering in a nuclear medium. The lowest-order 



contributions to the momentum imbalance are shown schematically in Fig. 14. At the lowest 
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Figure 13: Sketch for the di-jet mo- 
mentum imbalance measured at Fermilab 
E683. 



Figure 14: Sketch for the lowest order con- 
tributions to the averaged di-jet momentum 
imbalance in photoproduction. 



order, the di-jet momentum imbalance is approximately equal to the momentum imbalance 
between the two scattered partons which fragment into the two leading jets: ~ k^. As 
shown in Figs. 13 and 14, the averaged momentum imbalance (k^) is not exactly the same 
as the (fcy^) used in the experiment. The relation, assuming rotational symmetry, is simply 



(we always use rms averages) 



(16) 



Letting da^^{xp,i,p^) denote the photon-parton Born cross section with parton flavor {i = 
q, q, g) for fixed leading jet momentum £, we derived the double scattering contribution to 
the momentum imbalance [0, 



{k^Ei 



dH 



)4/3 



a. 



f do"^* 



i=q.q 



-Cg J dxfg/A{x)Ei-^{xp,i,p^) 



(17) 



with Cp = 4/3 and Cq = 3. The average value of k^ for all events in a region R of 
dijet phase space is found from Eq. ( [T7| ) by integrating over that region and dividing by the 
corresponding un- weighted cross section 0]. If we assume that asiix) in Eq- (P^ is evaluated 
at a typical value of the momentum transfer in R, and is kept approximately constant within 
R, the averaged momentum imbalance at the lowest order can be expressed as 



{kUR))4/3 = X'A'/'r 



a. 



CFafiR,p,) + CGafiR,p,) 



a]^{R,p^) + a]^{R,p^) 
In this expression, the un-weighted cross sections are defined as 



:i8) 



a. 



7A 



{R,Py) 



dH 



r da'^^ 
J dx fi/A (x) E£-^{xp, e, p^) . 



R Ee{2TT) 



(19) 



By comparing Eq. ([T8|) to data ||T6|, we found ~ 0.05 — 0.1 GeV^ |^. This value may be 
used to predict anomalous A-enhancement for other processes. 
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Figure 15: Sketch for the high or- 
der contributions to di-jets momen- 
tum imbalance in photoproduction. 



Figure 16: Sketch for the factorization of high or- 
der parton evolutions to jets from the lowest order 
parton level double scattering. 



The lowest-order calculation given above neglects the evolution of the final state partons 
into jets. We have argued above, however, that the evolution of the jet takes place outside 
of the nucleus. It is interesting to see how this decoupling of jet evolution occurs in our 
calculation. 

In perturbation theory, higher-order corrections due to the evolution of the final state 
partons in Fig. O arise from diagrams of the kind shown in Fig. |15|. The shaded subdiagrams 



may be thought of as "jet functions" , which represents the evolution of the outgoing partons. 
The first thing to note is that leading corrections due to final-state partonic bremmstrahlung 
are from small angle emission, which stays within the respective jets. Such collinear emission 
does not change the imbalance, fc^, although any of these fragmentation products can, in 
principle, exchange soft gluons with the nucleus. As pointed out in Ref. p|, however, soft 
gluon multiple scattering cannot resolve the details of the jets. The coupling of the final- 
state jet to soft gluons is coherent, and equivalent to a single "eikonal" line at leading power, 
whose direction and color are defined by the parton that initiates the jet. Thus, soft nuclear 



rescatterings factor from the jets, as shown in Fig. |T^. The combination of two such eikonal 
couplings, weighted by is enough to generate a lowest-order contribution to the matrix 
elements of Eqs. (|T2|) or (|13|). At yet higher orders, but remaining at first nonleading power 
in Q, additional soft gluon corrections generate nonabelian phases [|12|, (ordered exponentials 
of the gauge field), which serve to make the higher-twist matrix elements gauge-invariant, 
but do not otherwise affect A-dependence. One way of looking at this is that the di-jet kx 
imbalance is relatively insensitive to energy loss. Thus, the leading dependence of Eq. (|I]) is 
stable all orders in perturbation theory. 



3.2 ^4- Enhancement in Other Processes 

In order to test the theory, we need to identify other physical observables that are sensitive 
to the same twist-4 parton-parton correlation functions Ti{x,A) with i = q,g. One such 
process is direct photon production at measured transverse momentum, whose very moderate 
A-dependence has been measured by the E706 experiment at Fermilab. In Ref. [Q, it was 
found that the value of above, which produces a relatively large enhancement in dijet 
momentum imbalance, due to final-state interactions, produces a quite small 74-enhancement 
in photoproduction, due to initial-state interactions, consistent with experiment. This may 
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Figure 17: Lowest order diagram that 
contributes to the Drell-Yan trans- 
verse momentum broadening in hadron- 
nucleus colhsions. 



Figure 18: Sample double scattering am- 
plitudes that contribute to the NLO A- 
enhancement . 



shed some light on the long-standing observation that (initial-state) transverse momentum 
effects in Drell-Yan cross sections, to which we now turn, are also surprisingly small |1^, 113, 



19| 



In Ref . [0 , Drell-Yan transverse momentum broadening was calculated at the lowest order 



in pQCD. By evaluating the lowest order diagram in Fig. plus corresponding interference 
diagrams, it was found that the Drell-Yan transverse momentum broadening in hadron- 
nucleus collisions can be expressed in terms of the same twist-4 quark-gluon correlation 
function Tg(x, A) 

2 ^ (^^\ EqelJdx'f,/,ix')T,iT/x',A)/x' 
WtM/s 3 j j2^el!dx'f,/,ix')f,iT/x',A)/x' ^ ' 

where J2q runs over all quark and antiquark flavors, Cg is the quark fractional charge, and r = 
Q"^/ s, in terms of the lepton-pair invariant mass Q and hadron-hadron center of mass energy, 
^/s. Adopting the model in Eq. ([T^), the lowest order Drell-Yan transverse momentum 
broadening in Eq. (pOD can be simplified as 

{QlU/^ = (^] AM^/^ (21) 



By comparing Eq. (|21|) to data from Fermilab E772 and CERN NAlO experiments [|^, |T! 
it was found ^ that Drell-Yan data favor a A value: A^y ~ 0-01 GeV^, which is considerably 
smaller than the A^ extracted from the dijet momentum imbalance in photoproduction. 

It could be that the source of this discrepancy is experimental. On the other hand, we 
should also consider possible differences between the roles of multiple scattering in dijet mo- 
mentum imbalance and Drell-Yan. For dijet momentum imbalance in photoproduction, soft 
gluon multiple scatterings are purely final-state interactions, while the Drell-Yan transverse 
momentum broadening at the lowest order are caused by purely initial-state interactions. For 
initial-state interactions, there is strong interference beyond the leading order. Consider the 



diagrams in Fig. 0, which contribute to the Drell-Yan transverse momentum broadening, 

{Qt)a = A {QI)a=i + (Qr)4/3 (22) 
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at the next-to-leading order. Actually, these two diagrams can contribute to both terms 
in Eq. (p2D. In order to extract the nuclear size-enhanced contribution to the broadening 
(the A'^^^ term), we study the behavior of these diagrams at the poles near a: = 0. For the 
diagram in Fig. |T8|a, there are two potential poles labeled by A and B. When is small, 
contributions from these two poles cancel. That is, the NLO contributions to the Drell-Yan 
transverse momentum spectrum from this diagram do not have the dominant 1/Qt behavior, 
and therefore, its contributions to the {Qt)a/?, is much suppressed. When kx is small, the 
diagram in Fig. [T8| b is already included in the lowest order contribution to the broadening, 
because the gluon radiation is a part of the parton distribution fq/h- On the other hand, 
when kx is large, the only possible pole at B requires a finite value of a;, which ruins the 
A-enhancement ||^. 

From the above example, it seems plausible to us that the interference of initial-state 
radiation and multiple scattering suppresses ^-enhancement in Drell-Yan cross sections. By 
contrast, the corresponding effects of final-state radiation may cancel for an inclusive jet cross 
section, because we do not observe the fcr of radiation within the high-p^ jets. For Drell-Yan, 
the two terms of Eq. (^) are both from initial-state interactions, while in the dijet case, 
the first is primarily from the initial state, while the second is final-state. In the Drell-Yan 
case, destructive interference between the two mechanisms is possible, for di-jets, it is not. 
On this basis, we may expect a stronger A-enhancement to the dijet momentum imbalance 
in photoproduction than to the transverse momentum in Drell-Yan. Clearly, further study 
of this difference between initial-state and final-state interactions and related questions is in 
order. 



4 Conclusions 

In conclusion, we have argued that the nuclear size (or A^/^-type) enhancement caused by 
multiple scattering can be consistently calculated in pQCD, in terms of generalized factor- 
ization theorems. By studying enhancement from the nuclear medium, we can learn about 
strong interaction physics at twist-4. The first power correction measures new matrix el- 
ements {q{FF)q) and {F{FF)F). These matrix elements provide new insights into the 
nonperturbative regime of QCD. 

As reviewed above, the initial applications of pQCD to the first power corrections are 
interesting and generally successful. Different observables, such as jet broadening in deeply 
inelastic scattering |^ and pion transverse momentum broadening in deeply inelastic scatter- 
ing 1^, have been proposed and evaluated. The same techniques have been recently applied 
to many new physical observables M, 
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